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Cation Dependence of the Electronic States in Molecular Triangular Lattice System 
)8'-Z[Pd(dmit)2]2: A First-Principles Study 
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The electronic structure of an isostructural series of molecular conductors, yS'-Z[Pd(dmit)2]2, is systematically studied 
by a first-principles method based on the density-functional theory. The calculated band structures are fitted to the tight- 
binding model based on Pd(dmit)2 dimers on the triangular lattice. We find a systematic variation in the anisotropy 
of the transfer integrals along the three directions of the triangular lattice taking diff'erent values. The transfer integral 
along the face-to-face stacking direction of Pd(dmit)2 dimers is always the largest. Around the quantum spin liquid, X = 
EtMesSb, the other two transfer integrals become comparable. We also report sensible diff'erences in the distribution of 
wavefunctions near the Fermi level between the two dmit ligands of the Pd(dmit)2 molecule. 
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Electronic and magnetic properties of a series of anion 
radical salts, X[Pd(dmit)2]2 have attracted much attention re- 
cently.^ Here, monovalent cations X take Et3;Me4_3;Z, where y 
= - 2, Et = C2H5 , Me = CH3 , Z = P, As, and Sb, and dmit de- 
notes l,3-dithiol-2-thione-4,5-dithiolate. Most of the salts are 
isostructural at room temperature (RT) with the space group 
C2/c (called the /^'-type structure).^ These salts are paramag- 
netic Mott insulators at RT and exhibit a rich variety of ground 
states at low temperatures (LT) depending on the cation X. 
Many of them show antiferromagnetic (AF) order^ as in typ- 
ical Mott insulators, while others show interesting ground 
states, such as the valence bond solid (VBS) state"^ observed 
in the EtMesP salt, the quantum spin-liquid (QSL) state in the 
EtMcsSb salt,^~^ and the nonmagnetic charge-ordered (CO) 
state in the Et2Me2Sb salt.^'^ It was also reported that the 
EtMcsP salt^^'^^ and the salts showing the AF order become 
superconducting with the application of a hydrostatic or uni- 
axial pressure of less than 1 GPa.^^ To elucidate the origin of 
such a variety of ground states, it is necessary to clarify the 
cation dependence of the electronic structure. 

In this Letter, we report a systematic study on the electronic 
structure of y5'-X[Pd(dmit)2]2 by first-principles density- 
functional theory (DFT) calculations. The diflTerences be- 
tween the present results and the results of previous stud- 
ies obtained using the tight-binding model with the Htickel 
method (Huckel+TB)^'^'^^'^^ and DFT calculations 15 
be discussed. In addition, we present a detailed analysis of 
the character of Kohn-Sham orbitals (wavefunctions) near the 
Fermi level. 

The first-principles DFT calculations were performed by 
all-electron full-potential linearized augmented plane wave 
(FLAPW) method, which is one of the most reliable 
and accurate DFT methods. The exchange-correlation func- 
tional was treated within the generalized gradient approxi- 
mation (GGA) of the Perdew-Burke-Ernzerhof (PEE) func- 
tional..^^ Non- spin-polarized calculations were performed us- 
ing experimental crystal structures without geometrical opti- 



mization. 

Here, we work on nine members of yS'-type salts, which 
have the same space group C2/c at RT. Note that as for the 
EtMesP salt, the structure of a minor phase with C2/c was 
used, whereas its main phase showing the VBS state has 
a diflTerent space group Pli/m.^ As shown in Fig. 1(a), j3'- 
X[Pd(dmit)2]2 has a layered structure along the c-axis; the 
anion and cation layers consist of Pd(dmit)2 molecules and 
Et3;Me4_3;Z, respectively.^^ There are two equivalent anion 
layers in the unit cell, where Pd(dmit)2 molecules [Fig. 1(b)] 
stack face-to-face along the [110] and [110] directions with a 
strong dimerization (solid crossing column structure ^^). The 
Pd(dmit)2 dimers form an anisotropic triangular arrangement 
in each anion layer as shown in Fig. 1(c). 

Figure 2 shows the calculated band structures and Fermi 
surfaces of yS'-X[Pd(dmit)2]2, where X = (a) Me4P and (b) 
EtMesSb. Consistent with the previous reports, the bands 
crossing the Fermi level are isolated from the other bands and 
are 1/2-filled.^^'^^ Previous studies suggest that these bands 
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Fig. 1. (Color online) (a) Crystal structure of /5'-Me4P[Pd(dmit)2]2 (C2/c) 
at room temperature.^ Dotted lines represent the conventional unit cell, (b) 
Molecular structure of Pd(dmit)2. (c) Longitudinal view along the c-axis 
showing the ab-p\me. The interdimer transfer integrals tB, ts and tr are 
shown. 
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Fig. 2. (Color online) Band structure and Fermi surface of (5'- 
X[Pd(dmit)2]2 where X = (a) Me4P and (b) EtMesSb calculated by the first- 
principles DFT method. Dashed curves show the fitted results of the single- 
band TB model. The origin of the vertical axis with a dashed line shows the 
Fermi level. The present primitive unit cells (ap and bp) and the labels of 
/:-points are set to be the same as in the previous study. 



mainly originate from the antibonding pair of the highest 
occupied molecular orbitals (HOMO) of the two Pd(dmit)2 
molecules forming a (^imer.^'^^'^^'^^'^^ We will discuss this 
aspect and the character of the Kohn-Sham orbitals later. 

Now, let us focus on the cation dependence. Figure 3(a) 
shows the cation dependence of the bandwidth W of the bands 
crossing the Fermi level. The horizontal axis refers to the type 
of salt, i.e., X; we plot the salts according to the literature 
based on the Htickel + TB scheme. ^^'^^ The corresponding 
experimental ground states vary with decreasing Neel tem- 
perature to the QSL and CO states. Note that the magnetic 
structure of the jS'-type EtMesP salt at LT is unknown; thus, 
its order is determined by W. We find that W monotonically 
decreases as the number of ethyl (Ft) groups y increases, or 
the size of central atom Z in terms of the covalent radius 
becomes larger, from P, As, to Sb. This tendency is consis- 
tent with the chemical pressure eff'ect; larger cations (negative 
pressure) increase the volume of the unit cell and then, reduc- 
ing W. It is discussed that W can be considered as a control 
parameter of the electron correlation if the eff'ective Coulomb 
parameter does not show a marked cation dependence. On 
the other hand, it is reported that the anisotropy of the elec- 
tronic structure is correlated to the observed ground states.^ 
In this study, a characteristic cation dependence is seen in 
the anisotropy of the electronic structure shown in Fig. 2; the 
Fermi surfaces are composed of two equivalent cylinders as- 
sociated with two anion layers. By changing the cation from 
Ft2Me2Sb (No. 9) to Me4P (No.l), the cylinders are more dis- 
torted and less overlapped, indicating the electronic structure 
becomes more anisotropic. 

To evaluate the anisotropy quantitatively, we performed a 
tight-binding (TB) fit of the DFT bands and obtained the TB 
model using a single orbital of each dimer^' with three inter- 
dimer transfer integrals, tB, ts, and tr, in the Pd(dmit)2 layer, 
as defined in Fig. 1(c). An interlayer transfer integral across 
the cation layer, t^, is also included. This simple TB model 
well reproduces the DFT bands in all the cases (only two cases 
are shown in Fig. 2). The obtained transfer integrals with dif- 
ferent cations are plotted in Fig. 3(b). We find the following 
aspects: (1) is the largest. As the cation changes from left 



to right, (2) ts and tB decrease. (3) On the other hand, tr in- 
creases. The anisotropy represented by these transfer integrals 
well explains the cation dependence of the Fermi surfaces. 
Aspect (2) is consistent with the chemical pressure eff'ect on 
W, whereas aspect (3) is peculiar. This may be explained by 
the unique distortions of the Pd(dmit)2 molecule reported re- 
cently, which result in the shorter intermolecular distances be- 
tween S atoms along the diagonal direction even though the 
volume increases.^' 

For the materials on the left-hand side showing AF order, 
ts is greater than tr. However, as a result of (2) and (3), the 
order of the two integrals is reversed in the two right most 
salts, which show the QSL or CO state. Since the experimen- 
tal structure at LT is available for the EtMesSb (4.5 K^^) and 
Me4P salts (8 K^^), we have also calculated their band struc- 
tures and fitted the DFT bands to the same TB model (both 
maintain C2/c). Compared with the results obtained using the 
RT structure, tB and ts increase, whereas tr is almost the same. 
Consequently, the orders of the two transfer parameters tr and 
ts do not change, but they become close in the case of the 
EtMesSb salt. Aspect (1) still holds even though the LT struc- 
ture is used. Recently, Nakamura et al. have performed a 
DFT study on the EtMesSb salt with the LT structure and 
evaluated the transfer integrals between the Wannier orbitals 
constructed from the Kohn-Sham orbitals of the bands cross- 
ing the Fermi level. Their calculated parameter values are in 
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Fig. 3. (Color online) Cation dependence of (a) bandwidth W, (b) inter- 
dimer transfer integrals, tg^ts, and tr and interlayer transfer integral ^3 shown 
by solid circles, solid triangles, and solid squares, and crosses, respectively. 
The corresponding open marks indicate the transfer integrals obtained using 
the low-temperature structures, (c) Open circles show E\, the energy differ- 
ence between the HOMO and LUMO states of an isolated neutral Pd(dmit)2 
dimer, and solid circles show E2, the energy difference between the HOMO 
and LUMO-i-1 of the dimer, as shown in Fig. 4(a). 
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fairly good agreement with the present values derived from 
our TB fits. 

Let us summarize the present results in Fig. 3(b), in rela- 
tion to the experimentally observed properties. First, tB, the 
transfer integral along the direction parallel to the face-to-face 
staking of Pd(dmit)2 molecules (along [110] and [110]) is the 
strongest. This implies a one-dimensional (ID) structure as a 
first approximation. The ID chains are connected by ts and 
tr. In the salts showing AF order, ts is larger than t/, however, 
this diff'erence decreases in the compounds with lower Neel 
temperatures. Then, the two parameters become comparable 
around the two salts showing the competing ground states; 
the EtMesSb salt shows the QSL state and the Et2Me2Sb salt 
shows the CO state. Such a tendency indicates a possible 
relevance of eff'ective 1/2-filled electronic models or Heisen- 
berg models where ID chains are coupled via two frustrated 
interactions on an anisotropic triangular lattice. ^^'^^ 

Next, the present results are compared with the previous 
results obtained by Hiickel-hTB calculations using the same 
experimental structures.^ Aspects (2) and (3) agree with the 
previoues TB parameters between the antibonding HOMO of 
the dimers, whereas the aspect (1) diff'ers. Since tB and ts are 
nearly the same in each salt, tr/t, where t = {tB+ts )/2, was pro- 
posed as a key parameter to characterize the electronic struc- 
ture. Actually, the Htickel-hTB calculations for the LT struc- 
ture of the EtMesSb salt show a slightly larger value of tB than 
ts:tB = 34.7, ts = 32.6, and tr = 26.6 meV. However, the dif- 
ference is much smaller than the present results derived from 
the DFT calculations. 

Recently, Scriven and Powell have also employed first- 
principles DFT calculations and performed a similar system- 
atic study of y5'-X[Pd(dmit)2]2 by introducing the same ef- 
fective dimer model. They present five salts, except X = 
Me4As, Me4Sb, EtMesP, and EtMes As from the nine salts we 
have calculated. However, some of their TB parameters are 
diff'erent from ours. Their band structure of the EtMesSb salt 
is in good agreement with the present one, except that their 
position of the Fermi level is diff'erent from that in Fig. 2(b) 
and that of a previous study. ^^'^^ It is unclear why such dis- 
crepancies appear between their DFT-derived TB parameters 
and ours as well as those in ref. 15. 

Finally, we discuss the character of Kohn-Sham orbitals 
near the Fermi level. It was suggested that the bands at the 
Fermi level and the bands located around -0.4 eV mainly 
come from the antibonding pair of HOMOs and the bond- 
ing pair of the lowest unoccupied molecular orbitals (LU- 
MOs) of a Pd(dmit)2 molecule forming the dimer, respec- 
tively.^' ^^'^^'^"^'^^ We have calculated electron densities from 
the Kohn-Sham orbitals of these bands and found small but 
noticeable diff'erences in the distribution of the orbitals be- 
tween the left and right dmit ligands of the Pd(dmit)2 unit in 
the crystal. This was not recognized in previous studies. This 
disproportionation within the molecule is observed also in the 
local density of states (LDOS) of the left and right dmit lig- 
ands, as shown in Fig. 4(b). The contributions from the right 
are larger than those from the left at the bands located around 
-0.4 eV, while those from the left are larger than those from 
the right at the bands crossing the Fermi level. In these calcu- 
lations for the solid form, the disproportionation of the total 
charge density is small, almost negligible, compared with that 
of the charge density based on Kohn-Sham orbitals in each 



energy level. 

A similar disproportionation is also seen in the correspond- 
ing molecular orbitals (MOs) of an isolated Pd(dmit)2 dimer, 
but does not occur in an isolated monomer. The MOs and 
energy levels are calculated for one dimer/monomer in a su- 
percell, whose atomic positions are extracted from the crys- 
tal structures of the nine salts. The cation dependence of the 
HOMO-LUMO gap of an isolated monomer is almost negli- 
gible (-0.61 eV within the GGA), but the energy diff'erences 
between the dimer MOs (£"1 and E2) show a systematic cation 
dependence [Figs. 3(c) and 4(a)]. This agrees well with the 
energy diff'erence between the two bands of the solids, lo- 
cated at the Fermi level and around -0.4 eV. The two MOs 
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Fig. 4. (Color online) (a) Calculated molecular orbital levels for an isolated 
monomer and its dimer using a supercell. The geometries are adapted from 
y5'-Et2Me2Sb[Pd(dmit)2]2 at room temperature, (b) Local density of states 
(LDOS) of yS'-Et2Me2Sb[Pd(dmit)2]2. Dashed curves indicate the LDOS 
from the left ligand, and solid curves, that from the right dmit ligand. A ver- 
tical dashed line represents the Fermi level. The LDOS of the left ligand 
includes partial DOS from Cl-3, S3, S4, S7, S8, and SIO, while that from 
the right includes from C4-6, SI, S2, S5, S6, and S9 atoms. [Fig. 1(b)] 
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Fig. 5. (Color online) Contour plots of molecular orbitals for an isolated 
dimer, [Pd(dmit)2]2~, calculated using the (a) HOMO and (b) singly oc- 
cupied molecular orbital (SOMO) within GGA-PBE, and the (c) HOMO 
and (d) SOMO within B3LYP. The geometries are adopted from j3'- 
Et2Me2Sb[Pd(dmit)2]2 at room temperature. 
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of the isolated dimer also have noticeable differences in their 
amplitudes on the right and left dmit ligands, implying that 
the HOMO of the Pd(dmit)2 monomer hybridizes with the 
LUMO of the other monomer within the dimer. This is possi- 
ble since the HOMO-LUMO gap of the Pd(dmit)2 monomer 
is relatively smalP^ and the symmetry with respect to the left 
and right dmit ligands (the mirrors of C2v^^) is broken upon 
dimer formation in the solids. This molecular structure with 
the unique distortion is stable even after we optimize the in- 
ternal coordinates using the DFT method. 

It is also noted that pure DFT methods, such as GGA, of- 
ten underestimate such charge disproportionation. In fact, as 
shown in Fig. 5, our calculated MOs of the isolated dimer 
by the B3LYP method^^ using the GaussianOS code^"^ show 
that the degree of disproportionation of the corresponding 
MOs is much larger than that by the GGA.^^ We expect that 
this degree of freedom within the molecule (coupled with 
the molecular distortion) will be confirmed by future experi- 
ments and may play important roles in the stability of various 
ground states, especially for the CO state. ^ A detailed investi- 
gation using the fragment decomposition scheme^^ is now in 
progress and will be reported elsewhere. 

In summary, we have performed a systematic study of 13'- 
X[Pd(dmit)2]2 by first-principles DFT calculations and clari- 
fied the cation (X) dependence of the electronic structure. The 
dispersion of the DFT bands at the Fermi level is well repro- 
duced by the eff'ective TB model of the Pd(dmit)2 dimer on 
the anisotropic triangular lattice. The transfer integrals show 
that the parameter along the face-to-face stacking direction 
of the Pd(dmit)2 dimer ([110] and [110]), tB, is always larger 
than the other two parameters. We have also found that the 
other two parameters, ts and tr, become comparable in the re- 
gion where competing ground states (AF, QSL, and CO) are 
found experimentally. Furthermore, by analyzing the LDOS 
and Kohn-Sham orbitals of the solids, and the MOs of the iso- 
lated dimer, we revealed noticeable diff'erences in charge dis- 
tribution between the right and left dmit ligands of a Pd(dmit)2 
unit. This intramolecular degree of freedom may play impor- 
tant roles that clarify the various physical properties observed 
in the yS'-X[Pd(dmit)2]2 series. 
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